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^ ■ We describe inverse scattering for the matrix Schrodinger opera- 

, tor with general selfadjoint boundary conditions at the origin using 

t^H I the Marchenko equation. Our approach allows the recovery of the po- 

T-H ' tential as well as the boundary conditions. It is easily specialised to 

. inverse scattering on star-shaped graphs with boundary conditions at 

^ ' the node. 

> ■ 

m ; 

^ ; 1 Introduction 

O I The monograph by Agranovich and Marchenko [Ij provides the description 

Q I of the inverse scattering problem for the matrix Schrodinger operator on the 

^ ' semi-axis with Dirichlet boundary conditions at the origin. This followed 

CLc pioneering work by Gelfand, Jost, Krein, Levitan, Marchenko and other au- 

r-| ■ thors on the inverse problem. For a more complete description of the history 

! of this problem see the series of papers by Faddeev [191 ED] • 

^ I In this paper we generalise the work of Agranovich and Marchenko. We 

^ I consider the matrix Schrodinger operator with general selfadjoint boundary 

conditions at the origin. As a consequence the boundary conditions, in the 
form of a unitary matrix U, appear as a crucial element in our discussion 
of the direct and inverse problems. In particular we see that the scattering 
data, due to the continuous part of the spectrum, is of the form S{k) — U 
where S{k) is the scattering matrix and U is a matrix directly related to the 
boundary conditions (in fact it is the asymptotic of the scattering matrix). 
From our discussion of the inverse problem it is easy to see that the boundary 
conditions U along with the potential matrix are recovered in the course of 
the solution. 

The consideration of general boundary conditions, instead of Dirichlet bound- 
ary conditions, increases the complexity of the problem from an analytic and 
algebraic perspective. From an analytic point of view the characterisation of 
the set of scattering data becomes far more difficult. We are far from finding 
an equivalence between the set of permissible matrix Schrodinger operators 
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and a set of scattering data (on the other hand a complete characterisation is 
given in [T] for Dirichlet boundary conditions). In this paper we have concen- 
trated mainly on the algebraic complications that arise in the consideration 
of general boundary conditions. 

There has been a lot of work on systems with singular or finite rank pertur- 
bations [3], m [3l] and in particular they play an important role in the study 
of operators on graphs [HI [T21 HSl [H [ISl [13 [221 [2S]- Such systems may 
be useful in modeling the behaviour of real physical problems, in particular 
nanoelectronic devices jTEl UHl [23, 1291 [30], 126]. In these studies the choice 
of finite rank perturbation, or boundary condition, is an important part of 
the modeling of the dynamical system. On the other hand, as far as this 
author is aware, no work has been done on the inverse problem where the 
selfadjoint boundary conditions play a central role and are recovered as part 
of the solution. It is for this reason we believe that this paper will be of 
interest. 

As mentioned above many of the papers using finite rank perturbations in 
applications consider operators on graphs, rather than, as in our case, matrix 
systems. However, it is easy to see that the matrix system is a generalisation 
of the star-shaped graph which has an important, even central, role in ap- 
plications (see [in] where there is a discussion of modeling general quantum 
networks using such a system). 

There are many other approaches to the inverse problem aside from the ap- 
proach we have adopted here. There is an abundance of inverse spectral 
methods [HI [321 and the boundary control method [SI [3 [HI H] is another 
very useful approach. The boundary control method is particularly powerful 
when applied to inverse problems on finite graphs as it allows the recovery 
of the operator on the graph as well as the geometry of the graph [9]! The 
matrix Schrodinger operator is easily seen to generalise the Schrodinger op- 
erator on a graph with n rays and in [^ [2H] the authors consider the inverse 
problem on graphs using similar methods to the method used here (we note 
however that there are some errors in [21] , specifically the inverse problem is 
not correctly described). 

The contents of this paper are based on the results of the PhD thesis [21]. 

2 The Matrix Schrodinger Operator with Self- 
adjoint Boundary Conditions 

Let us consider the Hilbert space L^(]R_|_; C") consisting of the set of functions 
from ]R_|_ = [0, cx)) to C" satisfying 
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with associated inner product 

(/,^) = / fix)* ■g{x)dx. 



oo 





Here | ■ | is the usual norm in C" and ■* denotes the complex conjugate 

transpose for vectors and matrices. 

We define the matrix Schrodinger operator 

where Q = Q* is a hermitian matrix which satisfies 

[l + t)\Q{t)\dt<oo. (1) 







The norm, \A\, of a matrix acting on C" is defined here as the maximal 
eigenvalue of the matrix. The domain of £o is chosen to be the set of smooth 
functions with the (closed) support a compact subset of (0, oo) 

Dom(£o) = C^{^+] C") . 

In this case Cq has deficiency indices (n, n) — we are in the limit point case — 
and the selfadjoint extensions of £o can be parameterised by an r?, x n unitary 
matrix U using Neumann extension theory Here we take a slightly dif- 
ferent but equivalent approach [23] whereby the selfadjoint extensions are 
described by the boundary conditions at the origin 



(2) 



for some unitary U which is fixed in the remainder of the discussion. The 
latin subscript denotes differentiation. 

In order to simplify the presentation we assume that there are no virtual 
levels — ie. zero is not an eigenvalue — of C 



3 Properties of the Solutions 

Here we consider solutions Y{x^ k) of the eigenvalue equation 

CY = eY (3) 

for the most part ignoring boundary conditions and square integrability. 
It is convenient to define the involution 

Y^{x,k) = Y*{x,k) 
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as it is then clear that the Wronskian 

of two solutions of at the same value of k are constant (note the Wronskian 
as a function of $ and is a matrix valued hermitian symplectic form [23] • 
We also note that if we swop the order \E''''} is not constant). 

We define the standard solutions: 

Theorem 3.1 There exist solutions, and $, of the matrix equation ^ 
subject to the condition ^\) which satisfy the following boundary conditions 
at the origin 

lim^^o Q{x,k) =1, lima-^o Qx{x, k) = 
and are entire functions in the variable k. 

The proof of this theorem follows a standard argument (see theorem 1.2.1 
of Agranovich and Marchenko [1] for $, the discussion for B is similar [23]). 
We define the Jost solution using the asymptotic boundary condition 

lim F{x, k) = e^'^^I . 



Theorem 3.2 The Jost solution F{x,k) and its derivative Fx{x,k) are an- 
alytic in the upper half-plane of the variable k and continuous there, up to 
and including the real axis. 

This too is a standard result (it is a slight extension of theorem 1.3.1 in [1], 
see [211 for details). As a notational convenience we define 

F±{x,k) = F{x,±k) . 



Theorem 3.3 The Jost solution F^{x,k) can be written in terms of the 
transformation operator 



F+{x,k) = €'"''1+ K{x,t)e"''dt (4) 

J X 

with kernel K{x, t) which is bounded, absolutely integrable with respect to 
its second argument and the derivative K^{x,t) is also absolutely integrable 
with respect to its second argument. The kernel and the potential matrix are 
related by 

_^ dK{x,x) ^Q^^y 
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Again the proof may be found in [T], see mainly theorem 1.3.1. 

Finally we would like to define the scattered waves. We first define the 

function H(a;, k) as the solution of with boundary conditions 



(6) 



Here U is the fixed unitary matrix which defines the selfadjoint boundary 
conditions; A and B are defined above as a convenient shorthand. It is easy 
to see that A*B = B* A from which we immediately get that A;) satisfies 
the selfadjoint boundary conditions Clearly we can write S in terms of 
the standard solutions 

S(a;, k) = Q{x, k)A + ^{x, k)B , 

from which we see that H is entire in k, or in terms of the Jost solutions 

A;) = k)M^{k) + F+{x, k)M^{k) , 

where are some functions of the spectral parameter. Taking the Wron- 
skian 



W{FlE} 



FIB - fl^A 



±2ikM. 
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we get an expression for in terms of the Jost solutions. Here we use 

the constancy of the Wronskian and the two expressions for S in turn. The 
functions 

F±(fc) = F±(0,fc), f±,xik) = F±40,k) 
are known as the Jost functions. This gives us 



M, = ±— 
^ 2ik 



FIB 



^±,x^ 



(7) 



Choosing appropriate U we can make S equal the standard solutions. This 
allows us to write the standard solutions in terms of the Jost solutions. Eval- 
uating at the origin we get the following identities for the Jost functions 



F_Fl 



t 



F .Ft 



^+,x^+,x 



2ikl 
2ikl 

0. 



:t 

F-Ft -f/fV 
Note, these are not Wronskian relations. 

Lemma 3.1 The matrix M_{k) is invertable for finite real k. 
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Proof: We claim that for real k 



M^L* + LM* = I 

where L = FL ^.5 + FIA. Since k is real M*(/c) = m1(A;) and similarly 
L* = B*f_^x + A*F_ so the left hand side becomes 

[.4T_,^ - fi*F_] . 

Expanding this out and using AA* + BB* = I, AB* = BA* and the Wron- 
skian relations we immediately get the claimed equality. 
Suppose there is a non-zero a G ker (^M_ (j'-^^ for some A; G M. We have 
already shown that 

lini [a*M_{k)L*{k)a + a*L{k)M*_{k)a] =a*a^O, 

k^k 

which can only hold if L{k) has a pole at k. But this supplies a contradiction 
since, by theorem 13.21 the elements of the matrix L are bounded continuous 
functions for real finite values of k. □ 



2ik 



FIB - Fl „.A 



[ST_,, + A*F_] + 



Fl^B + FlA 



This generalises lemmata 2.2.2 and 2.4.1 of [T] where the result is proved 
for Dirichlet boundary conditions, ie —2kM_^ = ±Fj_. Consequently, for real 
k we define the scattered wave '${x,k) and scattering matrix S{k) in terms 
of k) 

^{x, k) = k)Mz\k) = F_{x, k) + F^{x, k)S{k) 

where 

(12) 

It is clear from the definition that the scattering matrix can be extended off 
the real axis as a meromorphic function and we show that away from any 
poles 

= s-^ . 

Taking the Wronskian of H"'' and S 

W{E\E} = A'B ~ B*A = , 

we see that it is zero. Furthermore, where the inverse Mr^ exists we can 
write S in terms of the scattering wave 

W{E\E} = MlW{^\<i/}M_ 

= 2ikMl [-1 + S^S] M_ = 

to get our result. In particular this gives the unitarity of the scattering 
matrix on the real axis. The high energy asymptotic of the scattering matrix 
is related to the selfadjoint boundary conditions by the following lemma: 



S{k) 



FlB 



FlB-Fl^A 
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Lemma 3.2 Given the self-adjoint C with associated unitary matrix U defin- 
ing the boundary conditions the scattering matrix has the asymptotic 



lim S{k) = U 



where U is a unitary hermitian matrix derived from U by applying the map 



to the spectrum of U. 

Here T is the unit circle in C. The proof follows from diagonalising U and 
the asymptotics of the Jost solutions. 

Our assumption that there are no virtual levels gives a simple proof of the 
fact that there are a finite number of discrete eigenvalues (see [1] pg. 38) 
although this statement is still true when there are virtual levels (as may 
be proved using Glazman's method of splitting, see [IJ theorem 2.1.1 for 
the Dirichlet case). Consequently, we have a finite number of negative (this 
follows from the existence of an integral equation representation of the Jost 
solutions which appears in the proof of therorem l3.2l) discrete eigenvalues kf. 
We write ki = iKi and choose the root ki > 0. The condition for a discrete 
eigenvalue is that, for some vector a, the square integrable vector function 
F^{x, ki)a satisfies the self-adjoint boundary conditions ([2D 



This is equivalent to det(M^(A;;)) = and, since ki is purely imaginary and 
Ml{ki) = -M*_{ki) we can write this as det(M_(A;;)) = 0. This gives: 

Theorem 3.4 The discrete eigenvalues of the self-adjoint Schrddinger op- 
erator correspond to the zeroes of 



in the upper half-plane. 

Theorem 3.5 The poles of MZ^ik) in the half-plane ^{k) > are simple. 

We say that a matrix has a simple pole at ki if it can be expanded as a power 
series in k — ki with the lowest order term (A; — ki)^^N^^i. 

Proof: Let us consider the entire solution H of the eigenvalue equation ([3]) 
and the solution S"'' of the 'adjoint' equation. We differentiate this adjoint 
equation with respect to k, multiply on the right by H, and subtract from it 
dH]) premultiplied by E^^. This gives 




det(M_(A;)) 




7 



We integrate the space variable from x to N and put k = ki to get 



xk^ 



N 



N 



2ki / [Ea]*Eadt 



(13) 



where a G ker (M_(iK/)) is non zero. Using the fact that a* eliminates 



a eliminates M_ and the constancy of the Wronskian the left hand side 
simplifies to 



a* Ml 



^-,xk^- 



N 



Since Q{ki) > all of the terms in the bracket are exponentially decreasing 
as — > oo so the upper limit vanishes leaving 



.F. 



M, a = 2ki / [Ea]*Ea dt . 



We now expand out M_,_ in terms of Jost functions and use the identities 
PHH) and M_a = to get 



POO 

ia*MiM_ i^a = / [Ha]*Sa dt ^ . 
Jo 

Now it is well known that MZ^(k) has a simple pole at k = ki iff the relations 

M_{ki)a = 
M_{ki)b + M_i^{ki)a = 

for some b implies that a = 0. Premultiplying the second relation by a*Ml 
gives 

a*MlM_b + a''MlM_f^a = 

and it is easy to see that the first term is zero: we use (IHl fTT!) and the fact 
that a* eliminates M|.. But this implies that a*M'^_M_ = which can 
only be true if a = 0. This implies that the pole is simple. □ 



4 The Inverse Scattering Problem 

In this section we derive the inverse scattering problem as an integral equation 
problem — the Marchenko equation (a discussion using the Riemann Hilbert 
problem is given in j24j). Our solution is to some degree formal; to treat 
this problem in its entirety we should first give a complete description of 
the space in which the scattering data exists and then show that the inverse 
scattering problem has a solution for any element of this space. A complete 
description of the space of scattering data for the matrix Scrodinger operator 
with general boundary conditions does not (to the knowledge of the author) 
exist. Indeed this space would almost certainly depend on the boundary 
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condition U (or what seems likely U). 

Nevertheless, given the condition ([T]) as well as the assumption of no virtual 
levels it is possible to show that the (scattering) data S{k) — is the Fourier 
transform of a hermitian matrix with integrable entries. This important 
result follows from theorem E.0.3 of [21] . We then have from theorem 3.4.1 of 
[1] that the Marchenko equation has a unique solution, again with integrable 
entries. We simply assume here a very narrow class of operators (certainly 
narrower than in [T]) but a class for which it is easy to show that the inverse 
problem has a solution. 



4.1 The Marchenko Equation 

The novel feature of the inverse problem for general boundary conditions is 
that the boundary conditions appear in the scattering data and the inverse 
problem through the asymptotic U. 

We start the derivation of the Marchenko equation by substituting (jll) into 
^ - F+U 

POO 

-^-F+JJ = e-'''''I + e'''%S -U)+ / K{x,t)e-'Hdt 

J X 

POO 

+ / K{x,t)e'''\S -U)dt. 



As noted above the Fourier transform 

fOO 



Gc{y)^^ P {S{k)-lj) e^'ydk 
exists (in fact is hermitian and integrable) so 

— / ('^{x,k) - F4x,k)u] e'^y dk = Gc{x + y) + K{x,y) + 



+ / K{x,t)Gc{t + y)dt (14) 

J X 

where we will only consider x < y. The left hand side is exponentially 
decreasing in the upper half plane (since x < y) and therefore we can close 
the contour of integration there. Since, theorem 13. 5[ the scattered wave has 
only simple poles at the eigenvalues there appears a sum of residues on the 
left hand side 

— ('^{x,k)-F+{x,k)U]e'^ydk = — E{x,k)MZ^e'^ydk 
27r y_oo ^ / 271- 

N 

■Y,^{x,ki)N_,ie''^^y 



N 

I 
1=1 



where N^^i is the residue of at the eigenvalue k = ki. As F+ is analytic 
in the upper half plane it makes no contribution. 
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To simplify the form of the residues we first need two auxihary results. Defin- 
ing Pi to be the orthogonal projection onto ker M\_{ki) we have 

E{x, k)U%F+ - iF+,,){k)Pi = F+(x, k)Pi . (15) 

This just follows if we evaluate the left hand side at x = 0, and the x 
derivative of the left hand side at a; = and see that we get the right hand 
side and its x derivative at x = 0. The second result 

M_/7^(F+-2F+,J(A:OPz = (16) 

follows if we expand M_ and use the left hand side of f|T5|) . again evaluated 

at X = 0, to get the Wronskian VF{f1, F+} which is zero. 

Let us consider the power series of M|(A;;) = — M*(fc;) and its inverse 

Mlik) = Ml{k) + {k-k)Ml,{ki) + --- 

where the subscript k denotes differentiation with respect to k and the sub- 
script / is an index of the zeroes ki. These expansions give us the following 
relations 

Miik)Nli = NlMiki) = (17) 
Ml,{k)Nl, + Ml{k)Ol, = Nl,Ml,{k) + Ol^^Mlik) = I . (18) 

From (ITTll we have PiN'!^_ ^ = ^ or taking the complex conjugate transpose 

N^,iPi = N^j . (19) 

We will now show that the left hand side of (11 5p is 'close to' the residue, 
namely 

U*iF+ - if+,.)Pi = iN.^iAiPi (20) 
where Ai is the positive definite hermitian matrix 

Ai= / FlF^{t,h)dt. 
Jo 

To see this we take the left hand side of fpU|) . multiply on the left by the 
complex conjugate transpose of equation ([T^ . expand out and use (fTBl) . then 
Mlik) Pi = to get 



F - F^ F 



p 



Using the same reasoning as was used to derive (1131) we see that the term in 
brackets F^ ^F+.a; — ^^F+ = 2kiAi which gives (pUj) . 
Going back to (ITSl) we get 

F+{x, ki)Pi = I S(x, ki)N_^iAiPi = I S(x, ki)N_,iBi 



10 



where Bi = PiAiPi + Pi'^ and we have used (fT9|) . It is clear that Bi is positive 
definite so we can define 

Q = PiB;'/^ . 

Then Cf = PiBf^ from which we finally get the desired form for the residue 

F+{x,ki)Cf = iE{x,ki)N.,i. (21) 

The Ci are known as normalisation matrices as the columns of F+{x, ki)Ci 
form a complete set of normalised eigenfunctions (this is a simple consequence 
of the above, for details see [T]). From (jlj) and ( fT^ we immediately get the 
Marchenko equation 



Theorem 4.1 Given the scattering data 

{S{ky, Ki, Gi,l = l,...,N} 

where S{k) is a unitary matrix, Gi are non-negative hermitian matrices 
and Ki are positive real numbers we can recover the potential of the matrix 
Schrodinger operator from the solution of the Marchenko equation 

poo 

G{x + y) + K{x,y) + K{x,t)G{t + y)dt = x<y, (22) 

J X 

where 

G{t) = y Gfe-^'' + — / (Sik) - U)e"''dk . 
Here U = limfc_>oo S{k). 

Proof: The potential is recovered from the solution of the Marchenko equa- 
tion using 1^. □ 



Corollary 4.1 Given the scattering data and the solution K{x,t) of the 
Marchenko equation we can recover the selfadjoint boundary conditions of 
the matrix Schrodinger operator from 

U = (^lo- z^,.|o)(^|o+ z^'.lo)-' . 

Proof: This follows from (j6]). The scattered wave can be found from the 
scattering matrix and (jl]). □ 
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4.2 The Diagonal Potential 

We finish with some brief comments on the case where Q{x) is a real diagonal 
matrix. This is easily seen to be equivalent to the scattering problem for a 
star-shaped graph with a single node and n semi infinite rays. In this case 
we expect some simplification to occur: 

Proposition 4.1 In the case of a diagonal potential the following scattering 
data is sufficient to recover the potential 

{Ri{k); Ki, -fi^i, i = 1, . . . ,n I = 1, . . . ,N} 

where Ri{k) = Sii{k), known as the reflection coefficients, are the diagonal 
elements of the scattering matrix and ji^i = Cf^^ are the diagonal elements 
of the squares of the normalisation matrices, known as the normalisation 
constants. 

Proof: In the case of a diagonal potential the kernel of the transformation op- 
erator K{x,t) is, like the Jost solution F+, a diagonal matrix. Consequently 
the diagonal elements of the Marchenko equation fl22l) form n independent 
scalar Marchenko equations. It is easy to see that these scalar Marchenko 
equations can be solved using only the above scattering data. □ 

There is good reason to expect that we can do better than this, but only 
if we are a priori given the form of the selfadjoint boundary conditions at the 
origin. We discuss this point for so called 'flux conserved' boundary condi- 
tions in the publication [23] and show there that it is possible to recover the 
potential with only n — 1 reflection coefficients and normalisation constants. 
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